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Abstract 

We study lepton flavor number violating rare B decays, b —» slflf, in a seesaw model with low 
scale singlet Majorana neutrinos motivated by the resonant leptogenesis scenario. The branching 
ratios of inclusive decays b —» sl±k T with two almost degenerate singlet neutrinos at TeV scale are 
investigated in detail. We find that there exists a class of seesaw model in which the branching 
fractions of b —> sry and r —> yy can be as large as 10 _1 ° and 10” 9 within the reach of Super B 
factories, respectively, without being in conflict with neutrino mixings and mass squared difference 
of neutrinos from neutrino data, invisible decay width of Z and the present limit of Br(y —> ej). 
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I. INTRODUCTION 


Thanks to the worldwide collaborative endeavor for B factories and neutrino experiments, 
our understanding of the flavor mixing phenomena in quark sector as well as in lepton sector 
has been dramatically improved over the past few years. Even though discoveries of neutrino 
oscillations in solar, atmospheric and reactor experiments gave a robust evidence for the 
existence of non-zero neutrino masses, we do not yet understand mechanisms of how to 


trac 


ive 


generate the masses of neutrinos and why those masses are so small. The most at 
proposal to explain the smallness of neutrino masses is the seesaw mechanism Q, 
in which super-heavy singlet particles are introduced. One of the virtues of the seesaw 
mechanism is to provide us with an elegant way to achieve the observed baryon asymmetry 


in our universe via the related leptogenesis js]]. However, the typical seesaw scale is of order 
10 10 ~ 10 14 GeV, which makes it impossible to probe the seesaw mechanism at collider 
experiments in a foreseeable future. Moreover, the leptogenesis at such a high energy scale 
meets a serious problem, the so-called gravitino problem, when it realizes in supersymmetric 
extensions of the standard model. Thus, it may be quite desirable to achieve the seesaw 
mechanism as well as the leptogenesis at a rather low energy scale. In these regards, scenarios 


of the resonant leptogenesis, in which singlet neu 
introduced, have been recently proposed 





, m 11 


h masses of order 1 ~ 10 TeV are 


, 12|, |13( Interestingly enough, the 


amplitudes of some flavor violating processes, which are highly suppressed in usual seesaw 
models with super-heavy singlet neutrinos, may be enhanced with such low scale singlet 
neutrinos. Thus, it is worthwhile to examine how we can probe a seesaw model with low 
scale singlet neutrinos at collider experiments and to find some experimental evidence for 
the seesaw model via probe of lepton flavor violating processes. 


Hi 


sl h l f 


In this paper, within the context of a low scale seesaw model based on SU (2) L xU(l) 
we study quark and lepton flavor violating (QLFV) rare decay processes such as b 
and lh —■► h'y, where lh and li denote heavy and light charged leptons, respectively. After 
the quark flavor changing neutral currents (FCNC) b —> sl + l~ had been discovered in B 


factories 


0 . 


it has been naturally expected that the next generation experiments could 


3 next gei 

probe well the lepton flavor violating (LFV) processes |lL Il6l . ll7L Il8 ] , which may eventually 
uncover the mechanism for the generation of small neutrino masses and the leptogenesis as 
well. Therefore, the precise predictions for those processes based on well motivated scenarios 
are very useful to find out if such scenarios can describe the nature correctly or not. Here 
we focus on the seesaw model motivated by the resonant leptogenesis scenario with two 









almost degenerate singlet neutrinos at TeV scale. In the context of the seesaw mechanism, 
lowering the singlet mass scale leads to an undesirable enhancement of the light neutrino 
rac neutrino couplings are guaranteed to be naturally small. However, 
, despite of the low mass scale of singlet Majorana neutrinos we can 
obtain light neutrino mass spectrum consistent with the current neutrino data by tuning 
the phase of the Yukawa-Dirac mass terms so that the two degenerate singlets contribute 
to the low energy effective Majorana mass terms destructively and the lepton number is 
approximately conserved. Interestingly enough, in such a scenario the sizable LFV processes 
and the suppression of the lepton number violation required in the effective Majorana mass 
for light neutrinos may coexist. We will show that the amplitudes of quark FCNC and 
LFV processes can be enhanced by considering some specific structure of Yukawa-Dirac and 
singlet Majorana mass matrices which are required to achieve the resonant leptpogenesis. 
We will obtain rather stringent constraints on QLFV processes by taking the constraints 
arisen from the invisible decay of Z boson, neutrino mass-squared differences and lepton 
flavor mixings measured at neutrino experiments and the experimental constraints of LFV 
processes such as lh —■► h'y- 

The paper is organized as follows: In section 2, we present the lepton flavor mixings 
of seesaw models with arbitrary number of singlet Majorana neutrinos. The analytical 
expressions for the branching fractions of the QLFV rare B decays are presented and the 
model independent bound on the branching fractions is obtained. In section 3, we study 
the low mass scale scenarios for singlet neutrinos and build the models in which the quark 
FCNC and LFV processes are enhanced. We give the predictions for the QLFV processes 
by taking into account the various constraints. In section 4, we summarize the results. 


masses unless the D 
as shown in 19 


II. LEPTON FLAVOR MIXIG OF SEESAW MODEL AND QLFV RARE B DE¬ 
CAY 


The lepton favor mixings of seesaw model are described in detail in Refs. [ 21 , 221. Here 
we extend the model to the case with arbitrary number of singlet neutinos and introduce a 
convenient decomposition of the Yukawa-Dirac mass term which is useful to our study. Let 
us start with a seesaw model described by following Lagrangian, 


r = 

J~~"m 


-yTLiN R J - ylLM - -N Rk M k N Rk + h.c., 


( 1 ) 
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and the neutrino mass matrix 


M u = 


0 m D 
mj) M 


( 2 ) 


where M is a N x N real diagonal singlet Majorana neutrino mass matrix and isa Dirac 
Yukawa mass term. Then, (3 + N) x (3 + N) neutrino mass matrix can be diagonalized by 
the mixing matrix V as, 


AfDiag = V ] M U V* 


(3) 


If the seesaw condition, i.e. ^ < 1, is satisfied, (3 + N) x (3 + N) unitary matrix V can be 
approximately parameterized as, 


V = 


M 


V mo 

~1h m D V { 1 n - ^rn D ^m D ^y 


(4) 


where V satisfies unitarity to the order of and In denotes an N x N unit matrix. Here 
3x3 submatrix V is not exactly unitary and can be written H as 


K= ( 1 - mD i^) Vo ' 

where Vo is a unitary matrix which diagonalizes the effective Majorana mass term, 

1 ^ 


(5) 


m eff = -m D j^rn D , 


VfjrrieffVo = Diag[ni,n 2 ,n 3 ], 

where ra* are the masses of three light neutrinos. 

For convenience, we introduce the following parametrization for rri£> 0 , 

/ rriDi 0 0 0 \ 

0 m D 2 0 0 


( 6 ) 


m D — fm D i, m D 2, •niDNj — I Ux u 2 


U N 


0 

V 0 


0 .. 0 
0 0 m DN ) 


where N unit vectors are introduced as, 


Ui = 


m DI 

m DI 


(7) 


( 8 ) 
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with mol = |rriDi|- By introducing the parameters of mass dimension Xj = (/ = 1 ~ 

N ), m e ff can be written as 

N 

m eff = ( 9 ) 

1=1 

The charged current is associated with 3 x (3 + N) submatrix of V. The deviation from the 
unitarity of 3 x 3 matrix V is given as 

3 N X 

T, ViaVja = % _ -J-T-UuUjj, 

a= 1 1=1 1 

E v “ v i = (10) 

i=e,/i,r 1=1 


where U = (u 1; u 2 ,u N ). 

Now we study the QLFV rare B decay processes b —> slflf in the context of the seesaw 
model we consider. We denote 4, as a heavy lepton and li as a light lepton, and the possible 
combinations for (if, If) are ('r ± , p T ), (/x ± , e T ) and (r^ e T ). By separating the contributions 
from the light neutrinos with masses n a (a = 1 ~ 3) and ones from the heavy neutrinos with 
masses M A _ 3 ( A = 4, 5, ...3 + N), we can express the amplitude for b —> sl^lf as 


rp _ V^Gfoiqed 

1 - ■ o X 


4vrs^ 


(TZ77 m Lw&) (utsy^Lvi) ^ ^ 


i=u c t 


a =1 


(3+AT) 

A=4 


( 11 ) 


where u b ,u s ,Uh and denote the spinor of bottom quark, strange quark, heavy lepton and 

3 ht anti-lepton respectively and L = x* = -S-, y a = -^§-, y A = i? is a Inami-Lim 

function presented as, 

E(x,y) = -(1 + x x M - /(y) + (1 - Zfy/'M - y/'M, 

4 x — y 4 x — y 


with /(x) = 


X log[x] 


( 12 ) 


x — 1 

We may neglect the up quark loop contribution because of the smallness of V ub Vf s . By using 
the unitarity relations for leptonic and quark sectors, 

3 3+iV 

E ^ = - E 


a=l A =4 

V tb v* = -v cb v*, 


(13) 
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we can simplify Eq. m as follows; 


py c 3+N 

T = -IL-9™ (W s ^Lu b ) {uH-fLvi) V tb V t * s ]T V, a V; a (E(x u x c , y A - 3)) , (14) 

47F A=4 

In Eq. ED, we only keep the contributions of the top quark, charm quark and heavy neutri¬ 
nos, and 


E(x t , x c , y A _ 3 ) = E(x t , y A - 3) - E(x t , 0) - E(x c , y A - 3) + E(x c , 0) 

_ /_3_ (xj - 8 x t + 4) log(ay) _ (y A _ 3 - 8 y A _ 3 + 4) log(?M- 3 ) A 

XtVA ~ 3 \4(1 _ x t ){l - y A - 3 ) 4(x t — l) 2 (xt — 2 M- 3 ) 4(|/ a _ 3 - l) 2 (y A - 3 - x t ) ) ’ 

(15) 


which agrees with the results in Refs. (l4l.l25|. The branching fraction of b —> sl h l+ can be 
easily obtained as 

It/ T/ 1 9 / \ 2 3+./V 

|V«>Vts| ^\m b ) ( OIqed\ 


Br(fe —> slh h + ) = Br(b —» crye)- 


|v c6 | 2 p( 


Sy SQ I5>,v^( 

\07T'S vy / A=4 


•P, Xc, DA—3) | , 
(16) 


where the phase factor P is given by P(x) = 1—8x 2 +8x 6 — x 8 — 24a: 1 log(x) and rrih denotes the 
heavier lepton (l^) mass. For numerical computation, we take Br(5 —> cv e e) = 0.107, m b = 
4.75 (GeV), m c = 1.25 (GeV), ckqed = m r = 1.78 (GeV), m /t = 0.106 (GeV), m w = 
80.4 (GeV). The branching fractions are then given as 


Br (b—>sr~e + ) = 1.0 x 10~' \S(r, e)| 2 , 

Br(6 —» sr~ y + ) = 1.0 x 10 -7 \S(r, /i)| J , 

Br(6 —■> sy~e + ) = 2.8 x 10~ 7 |S'( / u, e)| 2 , (17) 


where S(h,l ) is the suppression factor defined by 

3+N 

S(h,l) = ^2v hA Vi A E(x t ,x c ,y A _ 3 ), 

A=4 


N 


Xt 


^2 YT n hi n *uE{x t , x c , y /)• 
7=1 1 1 


(18) 


As can be seen above, the QLFV processes depend on u hI Ujj and Mj in E. 

If is not very small, there will be a chance to detect the QLFV processes at B factories 
near future. First, we can roughly estimate the branching fractions of the QLFV processes. 
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We notice that a model independent constraint on ^2f =1 can be o 
decay width of Z and charged curents lepton universality test |26|, |27|, 
coupling is suppressed compared with standard model prediction as, 


9 


2 cos 6\y 


Z ba Z»v b7/J (1 lb) v a 


stained from invisible 
28B because Z —»■ u a i7 b 


(19) 


a, 6=1 


where Z ba is related to the violation of the unitarity of V, 

Zi« = V Ok = - (v2u) u Zj-(u'v 0 ), a . ( 20 ) 

i=efir 

In the model under consideration, the effective number of light neutrinos N v is given by 

3 N 


Nu=Y, \Zba\ 2 = 3-2j2 

<2, 6=1 


7=1 


Xj_ 

M r 


( 21 ) 


From the experimental result j2fll |. N u = 2.984 ± 0.008, we can obtain the following bound, 


N 

E 

1=1 


Xj_ 

Mr 


= 0.008 ± 0.004. 


( 22 ) 


If the bound is dominated by N degenerate singlet heavy neutrinos with X 1 ~ X 2 X N , 
~ ^=, which is achieved in the case of low scale singlet neutrinos. 

From the fact that the lepton flavor universality of charged current is generally violated 


in seesaw model, one can obtain experimental bounds on e* = Y^i =1 \ u u\ 2 j^ (* = e , hi r ) 
from the lepton universality test of charged current interactions. The deviation from the 
universality is expressed in terms of the flavor dependent coupling g, defined as 

3 


9i_ 

g 2 


E v, 3 = 1 - «■ 


(23) 


< 2=1 


Using the dehnition of g t above, for instance, the decay width of W into charged lepton 
and neutrino is given by, 


E r \w^u 


Vn = 


g 2 M , 


w 


mf 2 mf 


a= 1 


487T 


E iwi 2 d--^-) 2 (i + 

1/1 w 


< 2=1 


2M2/’ 


g 2 M w mf 2 m\ 
48vr M^ )[L+ 2M^ h 


(24) 


where we ignore the neutrino masses. The experimental bounds on — e e and e r — e e was 
obtained from the ratios of the branching fractions of W decays and also from the ratios of 
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FIG. 1: The experimentally allowed region of (e^ — e e , e r — e e ) is shown, r, p denote the constraints 
obtained from r p leptonic decays. W denotes the bound obtained from W —> Uv decays. The 
invisible decay width constraints Eq. m is also shown for e e = 0. The dotted line coresponds to 
the prediction of class B model (NH) case. (See section IV in text.) 

the branching fractions of r —> fi(e)vv and p —> evv decays. The constraints obtained from 
r and p leptonic decays (W leptonic decays) are summarized in Eq.(13) (Eq.(6)) of Q, 

e M -e e = 0.0002 ± 0.0042 (0.002 ± 0.022), 

e T — e e = -0.0008 ±0.0044 (-0.058 ±0.028). (25) 

The off-diagonal elements of the correlation matrix of — e e and e T — e e are 0.51 and 0.44, 
respectively j^. In Fig.([TJ), we show the constrains Ea. (l251) on the plane (e^ — e e ,e r — e e ) 
by taking into account the correlations. In the figure, we also show the constraint obtained 
from Eq. m for the case that e e is vanishing. The constraint of Eq. 422) is written with e* 
as 


— € e) + ( e r ~ e e ) — (0.008 ± 0.004) — 3e e 


(26) 


From FiglU we can see that the constraints obtained from r and /i leptonic decays are 
consistent with Z invisible decay width constraint within lcr CL under the assumption e e = 0 
while the constraints obtained from W decays are not consistent and e T — e e < 0 seems to be 
required in this case. We come back to the lepton universality constraints when we consider 
the specific structure of Yukawa-Dirac mass term in the following sections. 

The Inami-Lim function E(-^p-) is shown in Fig. [21 The typical values for the Inami-Lim 









FIG. 2: Inami-Lim fiunctions —E (thick solid line) and 10 x (dashed line) as functions of the 
singlet neutrino mass M/y (GeV). 


function are \E\ ~ 2.8 ~ 8.1 for M = 200 ~ 2000 (GeV). Finally, from the fact that the 
factor UhiU* u depends on the flavor structure of m Dl the following relation 


2\uhiuu\ < \uhi \ 2 + \uu | 2 < ^ \uu\ = 1, (27) 

z=e,/x,r 


leads us to the constraint, \uhiUu[ : < 0.25. Then, numerical value of S(h,l ) for the case 
with N degenerate Mi and X t is 

|S(fc,()l 2 < (nF) E 

where we denote Mi = M{i = 1 ~ N ) and X{ r^j Xii — 1 ~ N). The upper bound 
of the branching fractions for b —* slf lf is roughly predicted to be 10“ n ~ 10 -10 for 
200 (GeV) < M < 2000 (GeV). 


M 2 


0.25 ~ 2 x 10 


-4 


2 x 10 


-3 


(28) 


III. LEPTONIC FCNC AND QLFV RARE DECAYS WITH LOW MASS SCALE 
SINGLET MAJORANA NEUTRINOS 

We now predict the branching fractions of the QLFV processes more concretely. As 
previously discussed, the branching fractions can be enhanced for rather large values of 
Xj/Mj. The large values of Xj/Mj are realized for low scale of Mj, which may not generally 
be consistent with neutrino data. The large values of Xj/Mj can be consistent with neutrino 
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(ni,n 2 ,n 3 ) 

Vo = (vi, v 2 , v 3 ) 

flavor dependence of b —► sl h 

Normal 

(°a/ A ™ sol> y/ Am atm + A "4l) 



Class A 

(0,|V 12 |,V 3 ) 

(W x u^,ui,u 3 )p 

Vh2 V[2 

Class B 

(0,V 3 ,|V 12 |) 

(u* x u*,u 3 ,ui)p 

V h 3V^ 3 

Inverted 

( \/ Am ltm - V Am atm, °) 



Class A 

(|V 12 |,V 3 ,0) 

(ui,u 3 ,u* x u *)p 

VhlV^ 

Class B 

(X 3 ,|X 12 |,0) 

(u 3 ,ui,u x x u* 3 )p 

Vh2 V[2 


TABLE I: The assignment of mass spectrum and MNS matrix. 


data when the contributions of Xj to m e ff in Eq. © are cancelled. Such a cancellation 
can be achieved by taking two almost degenerate small M 1j2 and tuning the relative phase 
between u x and u 2 so that those two terms contribute to m e ff destructively while keeping 
X 3 so small that its contribution to m e // is suppressed. Thus, we need some specific flavor 
structure of Yukawa-Dirac mass term in order to obtain an enhancement of the branching 
fractions. Let us assume X\ ~ X 2 3> X 3 so as for X i2 = X\ — X 2 and X 3 to be of order 
the light neutrino mass squared differences yj or yj Am^ tm . The relative phase of 
Yukawa-Dirac mass term from two singlet neutrinos N 3 and N 2 is tuned as u 2 = iu x . Then, 
m e ff becomes 

m eff = ~ ( u i u i Xi 2 + u 3 uJX 3 ) . (29) 

We further assume the orthogonality of u x and u 3 , i.e. u x • u 3 = 0, so that u x and u 3 can 

be directly related to MNS matrix. Then, there exists a massless state due to the alignment 

of u x and u 2 , which is assigned to n x = 0 for normal hierarchy and n 3 = 0 for inverted 

hierarchy. The other two masses are given by X 3 or Y 12 |. 

In Table [fl we classify the assignment of mass spectrum (n x , n 2 , n 3 ) and the unitary part 

of the mixing matrix Vo, where p is a diagonal Majorana phase which is irrelevant to the 

QFLV and LFV processes and thus we omit it from now on. Notice that Vo is identical 

2 

to the MNS matrix if we neglect its deviation from Wins, W — Wins = I 11 fact, 

2 

since the QLFV and LFV processes are already in the order of r - I ^ r at the leading order, the 
difference can be safely ignored. The flavor dependence of the amplitudes for the QFLV and 
LFV processes is then extracted in terms of the mixing angles of the neutrino oscillation 
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V /22V^ 2 

^13^sol^sol^atm 

0.33C13 

v r2V^ 2 

2 

^atm ^atm ^sol 

- 0.35 

v r2 v e2 

^13 SatmQsol ^sol 

— 0.33ci3 

v /23 v e3 

Cl3SatmSl3 exp(*5) 

0.71si 3 ci3 exp(zd) 

V r3V^3 

2 

Cl3<Sa,tmCatm 

0.5c? 3 

v r3 v e3 

Cl3C a tmSl3 exp (id) 

0.71si 3 ci 3 exp (id) 

V Ml V el 

C13 ^sol ^ sol ^at m 

—0.33ci3 

v rlV^ 

2 

^sol^atm^atm 

-0.16 

v rlV*! 

Cl3^sol^sol^atm 

0.33ci3 


TABLE II: The combinations of v,j relevant to the flavor dependence of QLFV and LFV decays. 


from Vo 


VmNS : 








/ 

C 13 CS 01 

£l3^sol 

0 N 


1 0 

0 

exp(— id) ^ 

^0 = 


'-’sol^atm 

Catm ^sol 

^13-5 atm 

+ Sl3 

-CsolSatm exp(i5) 

-SsolSat m exp(i(5) 

0 


V 

S sol $ atm 

^atm^sol 

^13^atm j 


^ -c so ic atm exp(i(5) 

-SsolCatm exp (i(5) 

° J 


(30) 


where we take s so i = 0.56, c so i = 0.84 and s atm = c atm = y 2 - 

In Table CD we present the relevant combinations of v,^ which correspond to the flavor 
dependence shown in the fourth column of Table I. The value of S 13 is very small and thus 
we ignore the terms of order 0(sl 3 ). Then, the suppression factor S(h,l) is approximately 
given by 


X - X - 

S{h,l ) ~ [—E(x t ,x c ,y 1 ) + —E(x t ,x c ,y 2 )\\i^l„ 


(31) 


where a denotes the index depending on the class and neutrino mass hierarchy, X\ ~ X 2 = 
X. And the term proportional to is not relevant at all and thus ignored. By using Eq. 
dm the ratios of the branching fractions are given by 

r±p=F'> * 2 


V V 


Br (6 —> 5 r =b e =F ) 

Br(& —■> st^ijiX) 

Br (6 —> s// ± e =F ) 

Br (6 -»• sr±/rF) P(^) 


p( > 
' mb > 


V M« V ea 


V V 


(32) 

(33) 


pt \ 
' td.l ' 


where P(^), and P(y^) are phase space factors and p^wy = 2.74 
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Class A NH (Class B IH) 


Class B NH 


Class A IH 


Br (b- 


r± P .T 


Br(6— >sr ± fi^) 


( ^atmCsol g sol j _ Q gQ 

Catm^atmC^ n j J 


Q Catm / 13 ) = 2.0 siq 
SatmCatm / - L ' 5 


^sol g sol g atm \ __ ^ ^ 

g snl ^atm^atm' , 


Br(b— ►s/z^e^) 
Br(6— >sr ± fi^ 1 ) 


^sol ^sol Catm \ 1 R __ 2 ^ 

Catm^atmC^ n ]' / ^V 


^atm^atn 


9; = 5.5*i 3 


c sol s sol c atm 

- - 

^snl Batin'S atm 


Q = 12.0 


TABLE III: Ratios of the branching fractions of b —> sl h 1+. P ^ = P(^) and P T = P(^)- 


Class A NH (Class B IH) 


Class B NH 


Class A IH 


Br(r— t-e'y) 
Br(r—>^7) 


^atm^atmC' 


sol 


■S at m C so l S sol -) =0.89 


Catmg13 ) = 2 . 0 sC 

S atm Catm / 1,3 


^sol^sol^atm 

~72 - 


sol 


Catm^atm 


= 4.5 


Br(^t—>e7) 
Br(r—>/i7) 


^sol ^sol ^atm 

- -5— 

Catm^atmC,^/ 


!_P_ 

Tr 


/ V 

y m r y 


= 5.0 


^atm^!3 

SatmCatm 


U± 

T t 


r£) =H4 


^-sol ^sol ^-atm 

- - 

g snl ^atm^atm 


LiL 

T t 


= 25 


TABLE IV: Ratios of the branching fractions of lh —> li y. 


In Table III, the numerical results of the ratios of the branching fractions given in Eqs. 
are presented. It can be seen that in Class B model for NH case, only b 


st /t 1 can 

be much larger than the other channel because of the absence of the suppression factor s 2 3 for 
r/i final states, while the branching fractions of the different channels in models except Class 


B for NH case are within a factor of 10. Furthermore, as discussed in Ref. jl4 1, there is strong 
correlation between QLFV processes and LFV radiative decays lh. —■► Z/y. Experimentally, 
there are stringent bounds as Br(/x —> ey) < 1.2 x 10 -11 [ 20 ], Br(r —> / 17 ) < 6.8 x ICC 8 
31] and Br(r -»■ erf) < 1.1 X 10- 7 (3.9 x 1(T 7 ) 32, 0 . The bounds on the LFV processes 
stringently constrain the branching fractions for QLFV processes. 

We note that the bound on -A from the invisible decay width of Z and the present upper 
bound on /i —> ey are not compatible with each other for Class A and Class B IH case, as 
shown in Table IV. The branching fraction for lh. —> Z/y is, 


Br(Z ft -> k y) = 


a 


ml m h 


|G| 2 . 


2567 r 2 s 4 w M 4 , T h . 

Numerically computing the pre-factors, the branching fractions are given by 


(34) 


Br(r —> /xy) = 5.4 x 1CT 4 |G(t,/ i)| 2 , 

Br(r —► ey) = 5.4 x 1CT 4 |G(t, e)| 2 , 

Br(yU —» ey) = 3.1 x 10“ 3 |G(/q e)| 2 , (35) 
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where G is a suppression factor defined by 


3 +N 


G = E Wafife„ 


a=4 

N 


^2 u u u h I jr F 2 (y I ), 


1=1 


Mr 


(36) 


_V+Vjs _ | as shown in 


where yj = and F 2 is Inami-Lim function, F 2 (y) = — Jq 
Fig. El 

With M x = M 2 = M < M 3 and X, ~ X, = X > X 3 , for Class A and Class B (IH) 
model, we predict, 

2 


Br(/i —> eq) = 3.1 x 10 


-3 


2X 

17 


(0.33) 2 F 2 (r/) 2 > 4.4 x 10 


-10 


(37) 


where we have used > 0.004 and M > 200 (GeV). Therefore, Class A and Class B (IH) 
models are excluded. If the bound on A- in Eq. (El is not taken into account, one can 
obtain from Eqs. ( 0351 ) . 


Br (b —> st~ /i + ) = 1.9 x 10 _4 Br(r —» /zq) 

Br(6 —■> sr~e + ) = 1.9 x 10 _4 Br(r —■> eq) 

Br (b —» s/i~e + ) = 9.0 x 10~ 5 Br(/i — > eq) 


irf * 13 x ir " 

|G(r/i )| 2 

{fl < 2.0 X 10-“ 
|G(re)| 2 

|S(,,e)|! < 1.1 x 10- 15 




G(Tfi) 

S(re) 


G(re) 
S(/j,e) 


(38) 


|G(/ze)| 2 G(/ie) 

While Eq. (EBt depends neither on the mass spectrum of heavy Majorana neutrinos nor on 
the flavor structure of Yukawa-Dirac mass terms, the ratio of ||)| depends on the details of 
them. For the present case with Mi = M 2 = M <C M 3 and X x = X 2 = X> X 3 , the ratio 
\S(h,l)/G(h,l )| is simply given as, 


S(h,l) 

2 

E(x t ,x c ,y) 2 

G(h,l) 


F 2 (y) 


= 98 (M = 200 GeV) - 285 (M = 2000 GeV). 


(39) 


Therefore, the upper bounds on the branching fractions given in Eq. m are translated to 


Br(6 —> st~ e + ) 
Br(6 —* st~ n + ) 
Br(6 —> sy~e + ) 


< (1.3 ~ 3.6) x 10" 9 

< (2.0 - 5.0) x 10~ 9 

< (1.1 - 3.0) x 10~ 13 


(40) 


The range of the upper bounds corresponds to M = 200 ~ 2000 (GeV). For Class A and 
Class B with IH case, by combining /i —> eq upper limit, we may obtain much tighter bound 
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on the other LFV and QLFV processes. Let us consider Class A (NH) and Class B (IH) 
cases. From Table ED and Table m the upper bounds on Br(/i —> ey) and Br(6 —> sfi e + ) 
obtained in Eq. m severely constrain the other modes, 


st e 


Br (6 
Br(6 —* st~ yU + ) 
Br(r -»• /zy) 
Br(r —> ey) 


< (0.41 ~ 1.1) x 1(T 13 , 

< (0.46 ~ 1.25) x 10“ 13 , 

< 2.4 x 10~ 12 , 


< 2.1 x 10 


-12 


(41) 


The bounds similar to Eq. HD are obtained for Class A IH case, 


Br(6 
Br (b 
Br(r 
Br(r 


sr~e + ) < (0.41 ~ 1.1) x 10- 13 , 
sr~n + ) < (0.92 ~ 2.5) x 10~ 14 , 
< 4.8 x 10“ 13 , 


m) 

■ ey) 


< 2.2 x 10 


-12 


(42) 


We note that the upper bound of the branching fractions of QLFV processes for Class A 
(NH, IH) and Class B (IH) models are 1CF 14 ~ 10 -13 and the branching fractions for LFV 
processes is 10” 13 ~ 10“ 12 . As we have already shown in Ea. (1.171) . the Class A model and 
Class B model for IH case can not satisfy the upper limit of the branching fraction of /i —» ey 
and constraint from the effective light neutrinos number N u simultaneously. This is because 
the former requires very small -A, while the latter requires larger jj. Below, we show Class 
B model for NH case may satisfy both constraints. Furthermore, the model predicts the 
large branching fractions of r —> /ry and b —> st/jl which are within the reach of near future 


Super B factories |3j,|35(. If we take into account the constraints coming from Br(/i —> ey) 
and the effective number of light neutrinos N„, we can have the parameter regions consistent 
with the present bounds only for Class B model with NH case. In this class, the stringent 
experimental limit on Br(/i —*■ ey) is not effective on Br(r —> /zy) and Br(6 — > st/i), because 
the former process is proportional to s 2 3 and thus ignorable, but the latter processes are not 
suppressed by the factor. 

In Fig. E we have shown the correlation of branching fractions between ji —» ey and the 
other LFV and QLFV processes in Class B model for NH case. The numerical results in 
Fig. E are obtained as follows: We first set X\ ~ X% = X and M 2 = M\ = M. From the 
constraint given in Eq. (1221) . the allowed range of M is , 


m D 1 


\[ 0 n 


V2 < M < 


m D 1 

\j ^ m i r 


V2, 


(43) 
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Br(i^ey) 


FIG. 3: Correlation between the branching fraction for fx —> ey and the branching fractions for 
b st fx (thick solid line) b —> sre (solid line) 6 — > s/ie (thin solid line), r — > /iy (long dashed line) 
r —>■ e 7 (dashed line) for Class B model with NH case. From left to right, the lines correspond 
to S13 = 0.02,0.05,0.1, repectively. The shaded region is excluded by the current bound on 
Br(/x —> ey). 

with a m in = 0.004, <r max = 0.012. When we fix mo 1, the allowed range of M is determined. 
By varying M within the above range, we plot the correlation between Br(/i —> ey) and the 
other five QLFV and LFV branching fractions. Here, S13 is a free parameter and is chosen to 
be 0.02, 0.05 and 0.1. mm is chosen to be 100 (GeV). As can be seen in Fig. |21 the present 
upper limit on Br(/i —»■ ey) gives a very tight bound on s 13 , typically smaller than 0.02. 
With this small s 13 , r —> ey and b —> sr(/i)e are also severely suppressed. Only b —> sr/i 
and r —> /ry are free from the suppression and the former branching fraction can be as large 
as 10~ 10 and the latter can be 10“ 9 . They are independent on small S13. 

We show in Fig. 0 and Fig. 0 the dependence of the branching fractions of b —> sr/i 
and r —/iy on mm and the heavy Majorana neutrino mass for the exact degenerate 
case, i.e. Mi = M 2 = M. We fix mm to 20,50,100,200 (GeV). Although the branching 
fractions become small as mm becomes small, the change of the branching fractions is 
within a factor 10. We also consider the non-degenerate case for Majorana neutrino masses 
(Mi eM / M 2 ) while keeping the degeneracy X\ ~ X 2 = X. By setting M 2 = RM\ and 
Mi = M, the dependence of the branching fraction b —> sr/i on the ratio R is studied. The 
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M GeV 


FIG. 4: Br (b —> st/i) vs. M for Class B (NH). From left to right, the curves correspond to 
uidi = 20,50,100,200 (GeV), respectively. 



0 1000 2000 3000 4000 5000 

M GeV 


FIG. 5: Br(r —> pi'y) vs. M for Class B (NH). From left to right, the curves correspond to 
rriDi = 20,50,100,200 (GeV),respectively. 

allowed range of the lightest heavy Majorana neutrino mass M of Eq. jSD is modified as 
+ T < M < ™ a D1 -y/l + Jy From Fig. |H1 we find that the lower and the upper 
limits of M become smaller as R becomes larger. However, the branching fraction does not 
change so significantly. 


IV. SUMMARY AND DISCUSSION 

As shown, the contributions of the singlet Majorana neutrinos to QLFV and LFV decays 
can be significant in the low scale seesaw model motivated by resonant leptogenesis. The 
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FIG. 6: Br(6 —> srfj.) vs. M for Class B (NH). We fix mpi = 100 (GeV). From right to left, the 
curves correspond to the ratio R = ^ 1, 2,10, respectively. 

branching fractions of inclusive decays b —> in the seesaw model considered in this 

paper depends on the suppression factor ^ which is arisen from the mixing between the 
singlet heavy neutrinos with three light neutrinos, and can be as large as about 10% without 
being in conflict with the neutrino mass squared differences from neutrino data and the 
current bound on the invisible decay width of Z boson. 

We have classified four classes of the model along with the light neutrino mass spectrum 
and the assignment of the mixing matrix Vo, and studied how the ratios of the branching 
fractions for the various channels of QLFV and LFV decays along with lepton flavors could 
be distinctively predicted in each class. We have found that only the class B for NH case 
presented in Table I survives the current limit on Br(fi —» ey) and the invisible decay 
width of Z boson. One may check if the model is consistent with the experiments of lepton 
universality tests. The class B for NH case predicts 

2X 

e e = —sl 3 < 5 x 10" 6 (si 3 < 0.02), 

2X 9 9 

% = ^-SatmOs - 0.004, 

2X 

€ r = -^-4m C L3 - 0-004. (44) 

The model predicts very small e e and = e T which are shown in Fig. © with dotted 
line. The model is consistent with the constraint of Z invisible decay width and the lepton 
universalty constraints from r and /j, decays while it may not be consistent with the lepton 
universality constraints determined by W decays. In this class, the branching fractions of 
b —* st /i and r —>• /ij are predicted to be as large as 10~ 10 and 10~ 9 , respectively. Such large 
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branching fractions can be tested in the future B factory experiments. The enhancment 
of the branching fractions of QFLV and LFV is originated from the one loop Feynman di¬ 
agrams in which the heavy Majorana neutrinos contribute to and it is non-susy contribution. 
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